\(\int \frac {x^3 \log (c (d+e x^2)^p)}{(f+g x^2)^2} \, dx\) [349]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (verified)
   Maple [C] (warning: unable to verify)
   Fricas [F]
   Sympy [F(-1)]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 25, antiderivative size = 155 \[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=-\frac {e f p \log \left (d+e x^2\right )}{2 g^2 (e f-d g)}+\frac {f \log \left (c \left (d+e x^2\right )^p\right )}{2 g^2 \left (f+g x^2\right )}+\frac {e f p \log \left (f+g x^2\right )}{2 g^2 (e f-d g)}+\frac {\log \left (c \left (d+e x^2\right )^p\right ) \log \left (\frac {e \left (f+g x^2\right )}{e f-d g}\right )}{2 g^2}+\frac {p \operatorname {PolyLog}\left (2,-\frac {g \left (d+e x^2\right )}{e f-d g}\right )}{2 g^2} \]

[Out]

-1/2*e*f*p*ln(e*x^2+d)/g^2/(-d*g+e*f)+1/2*f*ln(c*(e*x^2+d)^p)/g^2/(g*x^2+f)+1/2*e*f*p*ln(g*x^2+f)/g^2/(-d*g+e*
f)+1/2*ln(c*(e*x^2+d)^p)*ln(e*(g*x^2+f)/(-d*g+e*f))/g^2+1/2*p*polylog(2,-g*(e*x^2+d)/(-d*g+e*f))/g^2

Rubi [A] (verified)

Time = 0.16 (sec) , antiderivative size = 155, normalized size of antiderivative = 1.00, number of steps used = 10, number of rules used = 9, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.360, Rules used = {2525, 45, 2463, 2442, 36, 31, 2441, 2440, 2438} \[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=\frac {f \log \left (c \left (d+e x^2\right )^p\right )}{2 g^2 \left (f+g x^2\right )}+\frac {\log \left (c \left (d+e x^2\right )^p\right ) \log \left (\frac {e \left (f+g x^2\right )}{e f-d g}\right )}{2 g^2}+\frac {p \operatorname {PolyLog}\left (2,-\frac {g \left (e x^2+d\right )}{e f-d g}\right )}{2 g^2}-\frac {e f p \log \left (d+e x^2\right )}{2 g^2 (e f-d g)}+\frac {e f p \log \left (f+g x^2\right )}{2 g^2 (e f-d g)} \]

[In]

Int[(x^3*Log[c*(d + e*x^2)^p])/(f + g*x^2)^2,x]

[Out]

-1/2*(e*f*p*Log[d + e*x^2])/(g^2*(e*f - d*g)) + (f*Log[c*(d + e*x^2)^p])/(2*g^2*(f + g*x^2)) + (e*f*p*Log[f +
g*x^2])/(2*g^2*(e*f - d*g)) + (Log[c*(d + e*x^2)^p]*Log[(e*(f + g*x^2))/(e*f - d*g)])/(2*g^2) + (p*PolyLog[2,
-((g*(d + e*x^2))/(e*f - d*g))])/(2*g^2)

Rule 31

Int[((a_) + (b_.)*(x_))^(-1), x_Symbol] :> Simp[Log[RemoveContent[a + b*x, x]]/b, x] /; FreeQ[{a, b}, x]

Rule 36

Int[1/(((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))), x_Symbol] :> Dist[b/(b*c - a*d), Int[1/(a + b*x), x], x] -
Dist[d/(b*c - a*d), Int[1/(c + d*x), x], x] /; FreeQ[{a, b, c, d}, x] && NeQ[b*c - a*d, 0]

Rule 45

Int[((a_.) + (b_.)*(x_))^(m_.)*((c_.) + (d_.)*(x_))^(n_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*x)^m*(c + d
*x)^n, x], x] /; FreeQ[{a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && IGtQ[m, 0] && ( !IntegerQ[n] || (EqQ[c, 0]
&& LeQ[7*m + 4*n + 4, 0]) || LtQ[9*m + 5*(n + 1), 0] || GtQ[m + n + 2, 0])

Rule 2438

Int[Log[(c_.)*((d_) + (e_.)*(x_)^(n_.))]/(x_), x_Symbol] :> Simp[-PolyLog[2, (-c)*e*x^n]/n, x] /; FreeQ[{c, d,
 e, n}, x] && EqQ[c*d, 1]

Rule 2440

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))]*(b_.))/((f_.) + (g_.)*(x_)), x_Symbol] :> Dist[1/g, Subst[Int[(a +
 b*Log[1 + c*e*(x/g)])/x, x], x, f + g*x], x] /; FreeQ[{a, b, c, d, e, f, g}, x] && NeQ[e*f - d*g, 0] && EqQ[g
 + c*(e*f - d*g), 0]

Rule 2441

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))^(n_.)]*(b_.))/((f_.) + (g_.)*(x_)), x_Symbol] :> Simp[Log[e*((f + g
*x)/(e*f - d*g))]*((a + b*Log[c*(d + e*x)^n])/g), x] - Dist[b*e*(n/g), Int[Log[(e*(f + g*x))/(e*f - d*g)]/(d +
 e*x), x], x] /; FreeQ[{a, b, c, d, e, f, g, n}, x] && NeQ[e*f - d*g, 0]

Rule 2442

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))^(n_.)]*(b_.))*((f_.) + (g_.)*(x_))^(q_.), x_Symbol] :> Simp[(f + g*
x)^(q + 1)*((a + b*Log[c*(d + e*x)^n])/(g*(q + 1))), x] - Dist[b*e*(n/(g*(q + 1))), Int[(f + g*x)^(q + 1)/(d +
 e*x), x], x] /; FreeQ[{a, b, c, d, e, f, g, n, q}, x] && NeQ[e*f - d*g, 0] && NeQ[q, -1]

Rule 2463

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_))^(n_.)]*(b_.))^(p_.)*((h_.)*(x_))^(m_.)*((f_) + (g_.)*(x_)^(r_.))^(q
_.), x_Symbol] :> Int[ExpandIntegrand[(a + b*Log[c*(d + e*x)^n])^p, (h*x)^m*(f + g*x^r)^q, x], x] /; FreeQ[{a,
 b, c, d, e, f, g, h, m, n, p, q, r}, x] && IntegerQ[m] && IntegerQ[q]

Rule 2525

Int[((a_.) + Log[(c_.)*((d_) + (e_.)*(x_)^(n_))^(p_.)]*(b_.))^(q_.)*(x_)^(m_.)*((f_) + (g_.)*(x_)^(s_))^(r_.),
 x_Symbol] :> Dist[1/n, Subst[Int[x^(Simplify[(m + 1)/n] - 1)*(f + g*x^(s/n))^r*(a + b*Log[c*(d + e*x)^p])^q,
x], x, x^n], x] /; FreeQ[{a, b, c, d, e, f, g, m, n, p, q, r, s}, x] && IntegerQ[r] && IntegerQ[s/n] && Intege
rQ[Simplify[(m + 1)/n]] && (GtQ[(m + 1)/n, 0] || IGtQ[q, 0])

Rubi steps \begin{align*} \text {integral}& = \frac {1}{2} \text {Subst}\left (\int \frac {x \log \left (c (d+e x)^p\right )}{(f+g x)^2} \, dx,x,x^2\right ) \\ & = \frac {1}{2} \text {Subst}\left (\int \left (-\frac {f \log \left (c (d+e x)^p\right )}{g (f+g x)^2}+\frac {\log \left (c (d+e x)^p\right )}{g (f+g x)}\right ) \, dx,x,x^2\right ) \\ & = \frac {\text {Subst}\left (\int \frac {\log \left (c (d+e x)^p\right )}{f+g x} \, dx,x,x^2\right )}{2 g}-\frac {f \text {Subst}\left (\int \frac {\log \left (c (d+e x)^p\right )}{(f+g x)^2} \, dx,x,x^2\right )}{2 g} \\ & = \frac {f \log \left (c \left (d+e x^2\right )^p\right )}{2 g^2 \left (f+g x^2\right )}+\frac {\log \left (c \left (d+e x^2\right )^p\right ) \log \left (\frac {e \left (f+g x^2\right )}{e f-d g}\right )}{2 g^2}-\frac {(e p) \text {Subst}\left (\int \frac {\log \left (\frac {e (f+g x)}{e f-d g}\right )}{d+e x} \, dx,x,x^2\right )}{2 g^2}-\frac {(e f p) \text {Subst}\left (\int \frac {1}{(d+e x) (f+g x)} \, dx,x,x^2\right )}{2 g^2} \\ & = \frac {f \log \left (c \left (d+e x^2\right )^p\right )}{2 g^2 \left (f+g x^2\right )}+\frac {\log \left (c \left (d+e x^2\right )^p\right ) \log \left (\frac {e \left (f+g x^2\right )}{e f-d g}\right )}{2 g^2}-\frac {p \text {Subst}\left (\int \frac {\log \left (1+\frac {g x}{e f-d g}\right )}{x} \, dx,x,d+e x^2\right )}{2 g^2}-\frac {\left (e^2 f p\right ) \text {Subst}\left (\int \frac {1}{d+e x} \, dx,x,x^2\right )}{2 g^2 (e f-d g)}+\frac {(e f p) \text {Subst}\left (\int \frac {1}{f+g x} \, dx,x,x^2\right )}{2 g (e f-d g)} \\ & = -\frac {e f p \log \left (d+e x^2\right )}{2 g^2 (e f-d g)}+\frac {f \log \left (c \left (d+e x^2\right )^p\right )}{2 g^2 \left (f+g x^2\right )}+\frac {e f p \log \left (f+g x^2\right )}{2 g^2 (e f-d g)}+\frac {\log \left (c \left (d+e x^2\right )^p\right ) \log \left (\frac {e \left (f+g x^2\right )}{e f-d g}\right )}{2 g^2}+\frac {p \text {Li}_2\left (-\frac {g \left (d+e x^2\right )}{e f-d g}\right )}{2 g^2} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.07 (sec) , antiderivative size = 131, normalized size of antiderivative = 0.85 \[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=\frac {\frac {e f p \log \left (d+e x^2\right )}{-e f+d g}+\frac {f \log \left (c \left (d+e x^2\right )^p\right )}{f+g x^2}+\frac {e f p \log \left (f+g x^2\right )}{e f-d g}+\log \left (c \left (d+e x^2\right )^p\right ) \log \left (\frac {e \left (f+g x^2\right )}{e f-d g}\right )+p \operatorname {PolyLog}\left (2,\frac {g \left (d+e x^2\right )}{-e f+d g}\right )}{2 g^2} \]

[In]

Integrate[(x^3*Log[c*(d + e*x^2)^p])/(f + g*x^2)^2,x]

[Out]

((e*f*p*Log[d + e*x^2])/(-(e*f) + d*g) + (f*Log[c*(d + e*x^2)^p])/(f + g*x^2) + (e*f*p*Log[f + g*x^2])/(e*f -
d*g) + Log[c*(d + e*x^2)^p]*Log[(e*(f + g*x^2))/(e*f - d*g)] + p*PolyLog[2, (g*(d + e*x^2))/(-(e*f) + d*g)])/(
2*g^2)

Maple [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 2.08 (sec) , antiderivative size = 381, normalized size of antiderivative = 2.46

method result size
parts \(\frac {f \ln \left (c \left (e \,x^{2}+d \right )^{p}\right )}{2 g^{2} \left (g \,x^{2}+f \right )}+\frac {\ln \left (c \left (e \,x^{2}+d \right )^{p}\right ) \ln \left (g \,x^{2}+f \right )}{2 g^{2}}-p e \left (\frac {\munderset {\underline {\hspace {1.25 ex}}\alpha =\operatorname {RootOf}\left (e \,\textit {\_Z}^{2}+d \right )}{\sum }\left (\ln \left (x -\underline {\hspace {1.25 ex}}\alpha \right ) \ln \left (g \,x^{2}+f \right )-\ln \left (x -\underline {\hspace {1.25 ex}}\alpha \right ) \left (\ln \left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )}\right )+\ln \left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )}\right )\right )-\operatorname {dilog}\left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )}\right )-\operatorname {dilog}\left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )}\right )\right )}{2 g^{2} e}+\frac {f \left (-\frac {\ln \left (e \,x^{2}+d \right )}{2 \left (d g -e f \right )}+\frac {\ln \left (g \,x^{2}+f \right )}{2 d g -2 e f}\right )}{g^{2}}\right )\) \(381\)
risch \(\frac {\ln \left (\left (e \,x^{2}+d \right )^{p}\right ) f}{2 g^{2} \left (g \,x^{2}+f \right )}+\frac {\ln \left (\left (e \,x^{2}+d \right )^{p}\right ) \ln \left (g \,x^{2}+f \right )}{2 g^{2}}-\frac {p \left (\munderset {\underline {\hspace {1.25 ex}}\alpha =\operatorname {RootOf}\left (e \,\textit {\_Z}^{2}+d \right )}{\sum }\left (\ln \left (x -\underline {\hspace {1.25 ex}}\alpha \right ) \ln \left (g \,x^{2}+f \right )-\ln \left (x -\underline {\hspace {1.25 ex}}\alpha \right ) \left (\ln \left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )}\right )+\ln \left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )}\right )\right )-\operatorname {dilog}\left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =1\right )}\right )-\operatorname {dilog}\left (\frac {\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )-x +\underline {\hspace {1.25 ex}}\alpha }{\operatorname {RootOf}\left (e \,\textit {\_Z}^{2} g +2 \underline {\hspace {1.25 ex}}\alpha \textit {\_Z} g e -d g +e f , \operatorname {index} =2\right )}\right )\right )\right )}{2 g^{2}}+\frac {p e f \ln \left (e \,x^{2}+d \right )}{2 g^{2} \left (d g -e f \right )}-\frac {p e f \ln \left (g \,x^{2}+f \right )}{2 g^{2} \left (d g -e f \right )}+\left (\frac {i \pi \,\operatorname {csgn}\left (i \left (e \,x^{2}+d \right )^{p}\right ) {\operatorname {csgn}\left (i c \left (e \,x^{2}+d \right )^{p}\right )}^{2}}{2}-\frac {i \pi \,\operatorname {csgn}\left (i \left (e \,x^{2}+d \right )^{p}\right ) \operatorname {csgn}\left (i c \left (e \,x^{2}+d \right )^{p}\right ) \operatorname {csgn}\left (i c \right )}{2}-\frac {i \pi {\operatorname {csgn}\left (i c \left (e \,x^{2}+d \right )^{p}\right )}^{3}}{2}+\frac {i \pi {\operatorname {csgn}\left (i c \left (e \,x^{2}+d \right )^{p}\right )}^{2} \operatorname {csgn}\left (i c \right )}{2}+\ln \left (c \right )\right ) \left (\frac {f}{2 g^{2} \left (g \,x^{2}+f \right )}+\frac {\ln \left (g \,x^{2}+f \right )}{2 g^{2}}\right )\) \(523\)

[In]

int(x^3*ln(c*(e*x^2+d)^p)/(g*x^2+f)^2,x,method=_RETURNVERBOSE)

[Out]

1/2*f*ln(c*(e*x^2+d)^p)/g^2/(g*x^2+f)+1/2*ln(c*(e*x^2+d)^p)/g^2*ln(g*x^2+f)-p*e*(1/2/g^2/e*sum(ln(x-_alpha)*ln
(g*x^2+f)-ln(x-_alpha)*(ln((RootOf(_Z^2*e*g+2*_Z*_alpha*e*g-d*g+e*f,index=1)-x+_alpha)/RootOf(_Z^2*e*g+2*_Z*_a
lpha*e*g-d*g+e*f,index=1))+ln((RootOf(_Z^2*e*g+2*_Z*_alpha*e*g-d*g+e*f,index=2)-x+_alpha)/RootOf(_Z^2*e*g+2*_Z
*_alpha*e*g-d*g+e*f,index=2)))-dilog((RootOf(_Z^2*e*g+2*_Z*_alpha*e*g-d*g+e*f,index=1)-x+_alpha)/RootOf(_Z^2*e
*g+2*_Z*_alpha*e*g-d*g+e*f,index=1))-dilog((RootOf(_Z^2*e*g+2*_Z*_alpha*e*g-d*g+e*f,index=2)-x+_alpha)/RootOf(
_Z^2*e*g+2*_Z*_alpha*e*g-d*g+e*f,index=2)),_alpha=RootOf(_Z^2*e+d))+f/g^2*(-1/2/(d*g-e*f)*ln(e*x^2+d)+1/2/(d*g
-e*f)*ln(g*x^2+f)))

Fricas [F]

\[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=\int { \frac {x^{3} \log \left ({\left (e x^{2} + d\right )}^{p} c\right )}{{\left (g x^{2} + f\right )}^{2}} \,d x } \]

[In]

integrate(x^3*log(c*(e*x^2+d)^p)/(g*x^2+f)^2,x, algorithm="fricas")

[Out]

integral(x^3*log((e*x^2 + d)^p*c)/(g^2*x^4 + 2*f*g*x^2 + f^2), x)

Sympy [F(-1)]

Timed out. \[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=\text {Timed out} \]

[In]

integrate(x**3*ln(c*(e*x**2+d)**p)/(g*x**2+f)**2,x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=\int { \frac {x^{3} \log \left ({\left (e x^{2} + d\right )}^{p} c\right )}{{\left (g x^{2} + f\right )}^{2}} \,d x } \]

[In]

integrate(x^3*log(c*(e*x^2+d)^p)/(g*x^2+f)^2,x, algorithm="maxima")

[Out]

integrate(x^3*log((e*x^2 + d)^p*c)/(g*x^2 + f)^2, x)

Giac [F]

\[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=\int { \frac {x^{3} \log \left ({\left (e x^{2} + d\right )}^{p} c\right )}{{\left (g x^{2} + f\right )}^{2}} \,d x } \]

[In]

integrate(x^3*log(c*(e*x^2+d)^p)/(g*x^2+f)^2,x, algorithm="giac")

[Out]

integrate(x^3*log((e*x^2 + d)^p*c)/(g*x^2 + f)^2, x)

Mupad [F(-1)]

Timed out. \[ \int \frac {x^3 \log \left (c \left (d+e x^2\right )^p\right )}{\left (f+g x^2\right )^2} \, dx=\int \frac {x^3\,\ln \left (c\,{\left (e\,x^2+d\right )}^p\right )}{{\left (g\,x^2+f\right )}^2} \,d x \]

[In]

int((x^3*log(c*(d + e*x^2)^p))/(f + g*x^2)^2,x)

[Out]

int((x^3*log(c*(d + e*x^2)^p))/(f + g*x^2)^2, x)